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KINETIC  EFFECTS  IN  THE  CHEMISTRY  OF  DIAMOND  CVD  SOURCE 
GASES  AND  IMPLICATIONS  FOR  DIAMOND  GROWTH 

I.  INTRODUCTION 

The  growth  of  diamond  using  chemical  vapor  deposition  has  been  demonstrated 
for  a  wide  variety  of  source  gases.  Although  methane- hydrogen  remains  the  most 
common  growth  mixture,  in  recent  years  attention  has  focused  on  other  systems 
containing  carbon,  hydrogen,  and  oxygen.  Bachmann  et  al.  [1]  analyzed  a  large 
number  of  CVD  experiments  and  argued  that  it  is  not  the  specific  source  gases,  but 
rather  the  overall  relative  concentrations  of  the  atomic  species,  that  is  critical  in 
determining  whether  diamond  growth  is  possible.  The  mixtures  which  resulted  in 
diamond  growth  generally  possessed  C-H-0  stoichiometries  which  fell  within  a  well- 
demarcated  region  of  the  C-H-0  ternary  diagram  as  shown  in  Fig.  1.  Defining  the 
following  variable 

=  [C]I{[C]  +  [O]) ,  (1) 

they  had  determined  the  lower  boundary  separating  the  diamond  growth  and  no 
growth  regimes  to  be  given  by  the  line  Xc-s,  =  0.48.  The  upper  boundary  separating 
the  diamond  growth  and  non-diamond  carbon  growth  regimes  is  described  by  an 
arc  extending  from  nearly  the  CO  point  to  slightly  above  the  hydrogen  corner  of  the 
diagram.  Bachmann  et  al.  [1]  were  able  to  grow  diamond  from  CO  only  if  a  minimum 
of  0.5%  hydrogen  was  added  to  the  gas  mixture,  while  Saito  et  al.  [2]  could  produce 
only  amorphous  carbon  from  pure  CO.  The  shape  of  the  diamond  growth  region 
has  also  been  closely  approximated  by  Prijaya  et  al.  [3]  using  constrained  thermal 
equilibrium  calculations. 

These  results  can  be  understood  by  assuming  that  the  source  gases  have  reacted 
for  a  sufficiently  long  time  that  they  have  lost  essentially  all  knowledge  of  their  initial 
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state  and  that  a  quasi-thermodynamic  equilibrium  has  been  established.  The  orig¬ 
inal,  more  complex  molecules  have  been  reduced  primarily  to  a  mixture  of  atomic 
and  molecular  hydrogen,  small  hydrocarbon  species  such  as  CHn  and  C2Hm,  carbon 
monoxide,  carbon  dioxide,  oxygen  atoms,  and  hydroxyl  radicals.  For  gas  mixtures 
containing  equal  amounts  of  carbon  and  oxygen  in  thermodynamic  equilibrium,  the 
majority  of  the  carbon  flowed  into  the  reactor  is  tied  up  in  the  form  of  relatively 
non-reactive  carbon  monoxide.  The  stability  of  carbon  monoxide  is  so  high  that  only 
about  1%  is  dissociated  under  conditions  of  0.001  atm  and  4000  K,  for  example.  As 
the  mixture  is  made  more  oxygen  rich,  atomic  oxygen  and  hydroxyl  radicals  become 
prevalent  and  the  system  is  eventually  pushed  into  a  regime  in  which  not  only  is  dia¬ 
mond  growth  absent,  but  existing  diamond  [4]  and  non-diamond  carbon  [5]  are  etched. 
The  transition  from  diamond  growth  to  non-diamond  growth  on  the  carbon-rich  side 
of  the  growth  region  is  consistent  with  Goodwin’s  analysis  of  a  large  number  of  CVD 
experiments  [6].  He  observed  that  the  ratio  of  hydrogen  atom  flux  to  rate  of  carbon 
addition  must  exceed  a  critical  value  of  approximately  1  x  104-3  x  104  to  achieve 
diamond  growth.  Hydrogen  atoms  appear  to  be  important  both  in  creating,  through 
hydrogen  abstraction,  surface  radical  sites  where  growth  species  can  be  added  [7,8], 
and  in  “annealing”  the  surface  by  promoting  sp3  bonding  among  carbon  atoms  newly 
added  to  the  diamond  lattice  [8,9].  With  the  exception  of  the  very  earliest  work  in 
diamond  CVD  [10, 11],  all  successful  growth  recipes  have  included  large  amounts  of 
hydrogen. 

While  the  assumption  of  quasi-thermodynamic  equilibrium  may  be  valid  in  typical 
reactors  operating  at  relatively  low  flow  rates  and  pressures  of  about  30  Torr  and 
higher,  it  may  no  longer  be  true  under  operating  conditions  where  the  residence  time 
of  the  source  gases  in  the  reactor  is  smaller  than  or  comparable  to  the  time  required 
to  approach  thermodynamic  equilibrium.  In  particular,  our  interest  in  these  regimes 
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is  motivated  by  the  work  of  Rudder  et  al.  [12]  who  demonstrated  diamond  growth 
using  an  ethanol- water  rf-plasma  discharge  at  a  pressure  of  1  Torr.  Assuming  that 
the  gas-phase  chemistry  is  dominated  entirely  by  two-body  reactions,  the  residence 
time  required  for  the  source-gas  mixture  to  approach  equilibrium  scales  inversely 
with  pressure.  In  addition,  in  order  to  attain  higher  growth  rates  in  the  low  pressure 
processes,  it  would  be  advantageous  to  increase  the  mass  flow  rate  into  the  reactor 
thereby  increasing  the  total  amount  of  carbon  that  can  be  delivered  to  the  surface. 
The  combined  effects  of  low  pressure  and  high  flow  rate  can  lead  to  the  situation 
described  above,  where  the  transport  time  scales  are  short  compared  to  the  chemical 
time  scales. 

Although  the  results  described  in  this  paper  are  for  an  operating  pressure  of  0.001 
atm,  the  general  trends  observed  also  exist  at  higher  pressures.  We  had  shown  that 
the  characteristic  chemical  timescale  for  the  formation  of  atomic  hydrogen,  carbon 
monoxide,  and  growth  species  from  an  ethanol- water  mixture  is  inversely  proportional 
to  pressure  for  P  <  0.1  atm  [13].  In  addition,  the  residence  time  within  the  reactor 
is  proportional  to  the  pressure  divided  by  the  mass  flow  rate  m.  As  a  consequence, 
the  species  profiles  at  higher  pressures  ( Ph )  can  be  related  to  those  at  lower  pressures 
through  an  appropriate  change  of  scale  in  the  mass  flow  rate,  m  ->  m(P/Ph)2.  Since 
the  kinetic  effects  are  most  critical  for  the  low-pressure  processes,  the  focus  of  this 
paper  will  be  limited  to  results  obtained  for  P  =  0.001  atm. 

The  remainder  of  this  paper  is  organized  as  follows.  In  the  following  section,  we 
give  a  brief  description  of  the  elementary  reaction  mechanism  and  the  approach  used 
to  carry  out  the  perfectly  stirred  reactor  calculations.  The  results  of  the  calculations 
are  presented  for  six  source-gas  recipes,  comprising  two  different  C-H-0  stoichiome¬ 
tries,  in  section  III.  An  analysis  and  explanation  of  the  results  is  presented  in  section 
IV.  Finally,  in  section  V  we  summarize  our  findings. 
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II.  MODEL  FORMULATION 


The  elementary  reaction  set  used  for  the  calculations  described  in  this  paper  is 
based  on  a  combination  of  the  acetylene  combustion  mechanism  developed  by  Miller 
and  Melius  [14]  and  the  ethanol  combustion  mechanism  of  Curran  et  al.  [15].  The 
acetylene  mechanism  has  been  used  in  a  number  of  earlier  diamond  CVD  simulations 
(e.g.  Meeks  et  al.  [16]).  A  brief  description  of  the  mechanism  and  its  construc¬ 
tion  is  presented  below,  but  the  reader  is  referred  elsewhere  [13]  for  a  more  detailed 
treatment. 

The  mechanism  contains  a  total  of  249  reversible  reactions  involving  54  species. 
All  forward  reaction  rates  are  expressed  in  the  Arrhenius  form 

k  =  ATBe~E/RT ,  (2) 

and  include  enhanced  third-body  and  pressure  fall-off  forms.  The  acetylene  mecha¬ 
nism  was  chosen  as  the  starting  point  for  our  reaction  set  since  it  has  been  extensively 
validated  for  combustion  at  high  temperatures,  T  >  3000  K,  and  includes  a  detailed 
treatment  of  the  chemistry  of  species  such  as  C,  CH,  and  C2  which  are  important 
under  the  range  of  conditions  of  interest  in  this  study.  A  subset  of  the  reactions  from 
the  ethanol  combustion  mechanism  was  appended  to  the  acetylene  set  to  handle  the 
chemistry  of  ethanol  and  important  intermediate  species  not  explicitly  treated  in  the 
acetylene  mechanism.  This  subset  includes  23  reactions  directly  involving  ethanol 
and  five  additional  reactions  describing  the  decomposition  of  the  species  PC2H4OH, 
sC2H4OH,  CH3CHO,  and  C2H5O.  The  ethanol  reaction  set  was  originally  developed 
to  model  combustion  at  significantly  lower  temperatures,  T  <  2000  K,  but  for  the 
purpose  of  this  study  we  have  assumed  that  the  reaction  rates  in  this  mechanism  can 
be  extended  to  higher  temperatures. 

All  of  the  calculations  were  performed  using  the  Sandia  CHEMKIN  program 
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Perfectly  Stirred  Reactor  (PSR)  [17, 18].  The  basic  assumption  of  PSR  is  that  the 
mixing  of  reactants  is  so  complete  that  the  conversion  of  reactants  to  products  is 
totally  determined  by  the  chemical  reaction  rates  rather  than  diffusion,  convection 
or  other  transport  processes.  Fresh  source  gases  are  continuously  injected  into  the 
reactor  where  they  mix  and  react  with  the  older  material.  This  is  in  contrast  to 
plug-flow  reactor  simulations  where  there  is  no  mixing  of  new  and  old  material  and 
the  time-evolution  of  an  isolated  volume  of  the  reactants  is  calculated.  The  equations 
used  in  PSR  are  species  conservation 

Th(Yk-Yk*)-ukWkV  =  0,  (3) 


and  energy  conservation 

*h't(Ykhu-Ythl)  +  Q  =  0,  (4) 

k—1 

where  m  is  the  mass  flow  rate  through  the  reactor,  V  is  the  reactor  volume,  and 
Q  is  the  reactor  heat  loss.  The  subscripted  variables  T*,,  W4,  hk,  and  u are  the 
mass  fraction,  molecular  weight,  specific  enthalpy,  and  production  rate,  respectively, 
of  species  k.  The  starred  variables  represent  inlet  conditions  and  the  summation  in 
the  energy  equation  is  over  the  number  of  species  K.  The  mass  flow  rate  is  related 
to  the  residence  time  r  in  the  reactor  by  the  equation 

r  =  (5) 

m 

where  p  is  the  mass  density  in  the  reactor.  The  steady-state  conservation  equations 
are  solved  using  a  hybrid  Newton/time-integration  method.  The  PSR  program  can  be 
used  in  one  of  two  modes.  Either  the  temperature  and  pressure  of  the  reactor  can  be 
specified,  in  which  case  the  K  species  conservation  equations  are  solved,  or  the  reactor 
heat  loss  and  pressure  are  specified  and  the  K  +  1  species  and  energy  equations  are 
solved.  In  all  of  our  simulations  we  assume  that  the  reactor  temperature  and  pressure 
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are  known,  and  use  PSR  to  calculate  the  species  mole  fractions  in  the  reactor  for  mass 
flow  rates  corresponding  to  residence  times  ranging  from  1x10  4  s  to  1  x  103  s.  Since 
by  default  PSR  uses  the  results  of  a  thermodynamic  equilibrium  calculation  as  a  first 
guess  for  the  outlet  gas  composition,  it  is  easy  to  verify  that  the  PSR  calculations 
agree  with  thermodynamics  in  the  limit  of  long  residence  times.  For  nearly  all  of  our 
calculations  the  thermodynamic  equilibrium  conditions  are  closely  approached  for  a 
residence  time  of  1  x  103  s. 

In  addition  to  the  mole  fractions  of  individual  species,  we  are  also  interested  in 
the  flow  rate- dependent  behavior  of  several  derived  quantities  related  to  the  distri¬ 
bution  of  carbon  and  hydrogen,  the  rate  of  diamond  deposition,  and  the  quality  of 
the  deposited  material.  For  the  purpose  of  diamond  growth,  we  are  concerned  not 
just  with  the  total  concentration  or  mole  fraction  of  the  elemental  species,  but  with 
the  fractions  of  the  elemental  species  that  are  available  to  actively  participate  in  the 
diamond  growth.  Although  there  is  disagreement  over  which  hydrocarbon  species  are 
the  dominant  growth  precursors,  most  studies  suggest  that  it  is  the  Ci  radical  species, 
particularly  methyl  (CH3).  Diamond  growth  mechanisms  based  on  acetylene  as  the 
active  hydrocarbon  species  have  been  proposed  [19—24],  but  numerical  simulations  by 
Coltrin  and  Dandy  [25]  suggest  that  the  rate  of  diamond  deposition  due  to  acetylene 
is  generally  negligible  compared  to  that  for  carbon  atoms  or  methyl  radicals.  Ex¬ 
periments  and  modeling  by  Harris  and  Martin  [26]  suggest  that  the  growth  rate  and 
diamond  quality  are  significantly  lower  when  C2H2  is  the  active  growth  species  as 
compared  to  methyl  radicals.  Similar  conclusions  were  reached  by  Martin  and  Hill 
[27,  28].  Because  of  the  uncertainty  in  the  role  of  C2  species  in  diamond  growth,  we 
have  not  included  them  as  possible  diamond  precursors.  CH4  was  not  included  as  a 
growth  species  because  an  analysis  of  the  experiments  of  Chauhan  ct  dl .  [29]  by  Har¬ 
ris  [30]  suggests  that  although  CH4  can  participate  directly  in  diamond  growth,  the 
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surface  reaction  probability  for  CH4  is  approximately  six  orders  of  magnitude  lower 
than  that  for  CH3.  We  have  also  made  the  assumption  that  all  of  the  Ci  radicals  will 
play  the  same  role  in  diamond  growth.  Under  typical  CVD  conditions  where  the  flux 
of  hydrogen  atoms  to  the  surface  is  large  compared  to  the  C4  radical  flux,  C,  CH,  and 
CH2  species  which  are  added  to  the  surface  should  rapidly  have  their  dangling  bonds 
terminated  by  hydrogen  atoms. 

The  fraction  of  carbon  atoms  in  the  system  that  are  in  the  form  of  Ci  radicals 
can  be  defined  as 


/(Cirad)  = 


£3m=0X(CHm) 


(6) 


ELiEmEpnX(CnHroOp)’ 
where  the  quantities  X(C„HmOp)  are  the  mole  fractions  of  the  individual  species. 
The  quantity  X(CH2)  accounts  for  both  the  singlet  and  triplet  state  of  the  methylene 
radical.  The  summation  over  m  in  the  numerator  restricts  the  contribution  to  radical 
species.  In  the  denominator,  the  maximum  value  of  n  is  limited  to  six  since  our  reac¬ 
tion  mechanism  only  contains  up  to  C6  chemistry.  The  limits  for  the  summations  over 
m  and  p  are  assumed  to  be  consistent  with  the  number  of  carbon  atoms.  Similarly, 
the  fractions  of  carbon  in  the  forms  of  CH4,  C2,  and  heavier  (>  C3)  species  are  given 

by 

vrr’n.l 

(7) 


and 


,,CH  ,  = 

1  ''  Ej-iEmSp"X(C.HmO,)’ 

,  _  S»£»2X(C,HmO,) 

1(21  X^,,S-Sp»X(C.HmO „)’ 


\  _  Era=3  Em  Ep  ^A^CnHmOp) 
I{-  3)  ~  En=l  Em  Ep  nX(CnHmOp) 


(8) 


(9) 


Another  important  quantity  is  the  fraction  of  carbon  tied  up  in  the  form  of  carbon 
monoxide  and  related  species.  Once  carbon  reacts  to  create  carbon  monoxide,  it  is 
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essentially  inert  and  it  is  expected  to  play  very  little  role  in  diamond  growth,  at  least  as 
compared  to  active  radical  species.  Under  the  conditions  of  interest  using  the  reaction 
mechanism  described  earlier,  HCO  and  C02  undergo  reactions  to  form  exclusively  CO. 
CH30  leads  directly  to  the  creation  of  CH20  which  in  turn  decomposes  to  form  HCO. 
CH2OH  typically  reacts  to  form  CH20,  but  at  very  slow  flow  rates,  the  creation  of 
methyl  radicals  is  observed.  This  path  is  unimportant  though  since  the  mole  fractions 


of  CH2OH  at  these  flow  rates  are  less  than  1CT10.  The  fraction  of  carbon  in  the  form 
of  inert  species  then  is  given  by 


/(CO)  = 


X(CO)  +  X(HCO)  +  X(C02)  +  X(CH20)  +  X(CH2OH)  +  X(CH30) 

ELi  Em  Ep  nX(C„HroOp) 


(10) 


Finally,  we  define  the  fraction  of  hydrogen  in  the  form  of  active  hydrogen  atoms  as 

*(H) 


/(H)  = 


(11) 


En=0  Em  Ep  mX(C„HmOp) 

The  quality  of  the  deposited  diamond  can  be  predicted  from  the  relative  fluxes 
of  hydrogen  atoms  and  growth  species  to  the  surface.  An  analysis  of  the  conditions 
which  lead  to  the  growth  of  high-quality  diamond  indicates  that  a  minimum  of  ~ 
1  x  104-3  x  104  hydrogen  atoms  must  strike  the  diamond  surface  for  every  carbon 
atom  that  is  incorporated  on  the  surface  [6] .  If  we  neglect  boundary  layer  effects  and 
assume  that  all  CHm  radical  species  have  the  same  probability  of  sticking  to  a  surface 
radical  site,  we  can  make  the  following  estimate  of  the  ratio  of  hydrogen  atom  flux  to 
rate  of  carbon  atom  incorporation, 

X(H)/y/m^ 


9  = 


_ ,  (12) 

(ELo^CH^/v^b^’ 

where  7  is  the  sticking  probability  and  7i  is  the  fraction  of  surface  sites  that  are 

open.  We  used  the  values  of  7  =  0.33  based  on  the  surface  chemistry  mechanism  of 
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Coltrin  and  Dandy  [25]  and  Z  —  0.085  derived  from  the  hydrogen  atom  abstraction 
and  recombination  probabilities  of  Brenner  et  al.  [31,32]. 

A  maximum  theoretical  rate  of  diamond  deposition  can  be  defined  as 

Dmax  =  myc/(Cirad),  (13) 

where  Yq  is  the  total  carbon  mass  fraction  of  the  inlet  gases.  The  dimensions  of  Dmax 
are  (mass/time)  and  were  chosen  in  this  way  rather  than  in  the  usual  dimensions 
of  (length/time)  so  that  the  results  will  be  independent  of  reactor  geometry  and 
growth  area.  Our  expression  for  Dmax  is  based  on  the  assumption  that  Ci  radicals 
are  converted  to  diamond  with  100%  efficiency.  Obviously,  this  is  an  unrealistic 
assumption  since  it  neglects  the  growth  of  non-diamond  carbon,  the  finite  transport 
rates  of  the  growth  species  and  hydrogen  atoms,  and  the  non-unity  probability  of 
incident  growth  species  sticking  to  the  surface.  Nonethless,  our  definition  is  still 
useful  for  several  reasons.  First,  it  sets  an  absolute  upper  bound  on  the  diamond 
deposition  rate.  Second,  taken  together  with  the  quality  factor  defined  above,  it 
can  be  used  to  help  establish  optimal  flow  rates  and  reactor  operating  conditions. 
Finally,  although  Dmax  may  be  substantially  larger  than  deposition  rates  that  can 
be  realistically  expected,  it  allows  us  to  make  at  least  rough  comparisons  between 
different  source-gas  recipes. 


III.  RESULTS 

In  this  section  we  present  the  results  of  our  calculations  for  a  variety  of  CVD 
source-gas  recipes.  Three  of  the  mixtures,  ethanol-water-hydrogen  21:20:1,  methane- 
oxygen  84:41,  and  acetylene-oxygen-hydrogen  42:41:126,  all  possess  the  same  overall 
C-H-0  stoichiometry  of  42:168:41.  These  mixtures,  which  will  be  referred  to  as  CO- 
rich  systems  since  they  all  contain  large  amounts  of  carbon  and  oxygen  in  nearly 
identical  proportions,  are  labeled  as  “mixl”  on  Fig.  1.  This  choice  of  stoichiometry 
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was  motivated  by  our  interest  in  the  ethanol-water  system  [12,  13,  33]  and  is  just 
slightly  carbon-rich  relative  to  pure  methanol  (C-H-0  =  1:4:1).  A  second  set  of  three 
mixtures,  ethanol-hydrogen  1:81,  methane-oxygen-hydrogen  4:1:160,  and  acetylene- 
oxygen-hydrogen  2:1:166,  all  have  an  identical  C-H-0  composition  of  2.168.1,  with 
the  stoichiometry  chosen  so  as  to  be  in  the  range  of  typical  methane-oxygen-hydrogen 
recipes.  This  second  set  of  mixtures,  labeled  as  "in ix2  on  Fig.  1,  will  be  referred 
to  as  CO-lean  systems  since  carbon  monoxide  is  a  relatively  minor  component  of  the 
fully  reacted  gases.  Our  choice  of  nomenclature  is  somewhat  arbitrary,  but  it  does 
provide  a  descriptive  way  of  differentiating  the  two  sets.  All  six  of  the  mixtures  have 
an  H/(C-0)  ratio  equal  to  168.  By  choosing  different  gas  mixtures  with  identical 
C-H-0  ratios,  we  are  able  to  distinguish  effects  due  to  the  specific  source  gases  from 
those  due  to  the  overall  stoichiometry.  The  H / (C-O)  ratio  was  chosen  to  be  a  constant 
in  order  to  confirm  that  the  quality  factors  for  all  six  cases  converge  to  the  same  value 
in  the  limit  of  low  flow  rates  (long  residence  time). 

The  carbon  and  hydrogen  distributions,  diamond  quality  factors,  and  maximum 
theoretical  deposition  rates  are  displayed  in  Figs.  2-7  for  the  six  recipes  described 
above.  In  all  of  the  simulations,  the  reactor  volume  is  1000  l,  the  temperature  is  3000 
K  and  the  pressure  is  0.001  atm.  For  other  reactor  volumes,  the  scaled  flow  rate  for 
a  1000  l  reactor  is  related  to  the  actual  flow  rate  by 


^1000  — 


1000 

V 


-m. 


(14) 


The  maximum  deposition  rates  for  two  reactors  operating  at  the  same  m/V  ratio  are 


related  by 


D 


max 


V 

1000 


D 


maxlOOO* 


(15) 


For  all  six  cases,  Dmax  drops  with  decreasing  flow  rate  for  low  flow  rates,  with 
the  dependence  becoming  linear  in  the  limit  rn  — >  0.  This  reflects  a  convergence  to 
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thermodynamic  equilibrium  in  the  reactor  vessel  as  r  — »  oo.  The  deposition  rate 
is  directly  proportional  to  the  flow  rate,  and  depends  only  on  the  overall  C-H-0 
stoichiometry.  At  high  flow  rates,  the  behavior  of  Dmax  depends  sensitively  on  the 
nature  of  the  carbon  source.  Feed  gases  containing  ethanol  or  acetylene  both  result 
in  a  well  defined  maximum  in  Dmax,  while  source  gases  relying  on  methane  as  the 
carbon  source  lead  to  a  deposition  rate  that  asymptotically  approaches  a  maximum 
in  the  limit  m  oo. 

The  behavior  of  the  quality  factor  is  slightly  more  complicated  since  it  depends  on 
both  /(Cirad)  and  /(H).  In  the  limit  m  — >  0,  the  quality  factors  for  the  six  cases  are 
independent  of  m  and  all  converge  to  the  same  value,  as  shown  in  Fig.  8.  \F  exhibits  a 
maximum  in  the  vicinity  of  m  =  10-2  g/s,  but  the  maxima  for  the  CO-lean  mixtures 
are  considerably  more  pronounced  than  those  for  the  CO-rich  mixtures.  The  most 
striking  differences  in  the  quality  factors  occur  at  high  flow  rates.  With  the  exception 
of  the  CO-lean  acetylene  mixture,  'F  drops  rapidly  as  the  flow  rate  is  increased.  The 
quality  factors  for  the  CO-lean  ethanol-hydrogen  and  acetylene  mixtures  show  local 
minima  at  high  flow  rates,  while  the  remaining  mixtures  have  quality  factors  that 
drop  monotonically  with  increasing  m.  The  details  of  this  behavior  are  probably 
unimportant  since  they  occur  in  flow  regimes  where  \F  is  well  below  that  required  for 
diamond  growth. 

The  carbon  and  hydrogen  distributions  show  that  in  the  limit  rh  —*  0,  the  inlet 
gases  are  converted  almost  entirely  to  carbon  monoxide,  atomic  hydrogen  and  Ci  rad¬ 
ical  species  as  predicted  by  thermodynamics.  At  the  other  extreme  of  high  flow  rates, 
the  inlet  gases  pass  through  the  reactor  with  minimal  conversion  to  the  thermody¬ 
namically  stable  products.  In  all  cases,  /(CO)  and  /(H)  increase  steadily  as  the  flow 
rate  is  decreased  and  /(>  C3)  exhibits  a  global  maximum  for  intermediate  flow  rates. 
The  behavior  of  /( C2)  and  /(CH4)  depend  on  the  form  of  the  carbon  source.  For 
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mixtures  containing  ethanol  or  acetylene,  /( C2)  decreases  steadily  as  m  — *  0  while 
/(CH4)  reaches  a  maximum  for  intermediate  flow  rates.  The  characteristic  behaviors 
of  /( C2)  and  /(CH4)  as  a  function  of  m  are  exchanged  when  methane  is  used  as  a 
carbon  source.  The  most  unusual  trend  is  that  shown  by  /( Cirad)  as  the  flow  rate 
is  varied.  At  first,  /(Cirad)  rises  steadily  as  the  flow  rate  is  decreased  from  infinity, 
but  it  passes  through  a  global  maximum  and  then  a  local  minimum  before  settling  to 
the  value  predicted  by  thermodynamic  equilibrium  at  low  flow  rates. 

IV.  ANALYSIS 

Based  solely  on  the  behavior  of  Dmax ,  the  strategy  for  diamond  growth  in  a 
low-pressure,  hot-filament  or  microwave-plasma  reactor  would  appear  to  be  straight 
forward:  for  source  gases  containing  methane,  increase  the  flow  rate  until  the  diamond 
deposition  rate  levels  off;  for  recipes  utilizing  ethanol  or  acetylene,  increase  the  flow 
rate  until  the  predicted  maximal  growth  rate  is  obtained.  Unfortunately,  two  factors 
related  directly  to  the  gas-phase  chemistry  severely  limit  the  potential  growth  rates. 
First,  the  quality  of  the  deposited  material  depends  on  the  ratio  of  atomic  hydrogen 
to  Ci  radicals.  In  most  of  the  cases  examined,  with  the  notable  exception  of  the 
CO-lean  acetylene  mixture,  the  Ci  radicals  are  created  more  readily  than  hydrogen 
atoms  at  high  flow  rates,  leading  to  unfavorable  values  of  T.  Second,  it  is  believed 
that  heavier  hydrocarbon  species  may  be  responsible  for  the  growth  of  non-diamond 
carbon.  While  there  is  uncertainty  over  the  role  of  C2  species  in  diamond  and  non¬ 
diamond  carbon  growth,  it  is  likely  that  the  heavier  hydrocarbon  species  (>  C3)  do 
not  play  a  role  in  diamond  growth  and  may  in  fact  lead  to  the  deposition  of  non¬ 
diamond  carbon.  Particularly  for  the  CO-rich  mixtures,  the  fraction  of  carbon  in  the 
form  of  C3  and  heavier  species  can  be  significant  (~  10%)  and  considerably  greater 
than  the  fraction  in  the  form  of  Ci  radicals. 
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While  the  hydrogen  and  carbon  distributions  in  the  limit  of  slow  flow  rates  are 
specified  entirely  by  thermodynamics,  the  distributions  at  higher  flow  rates  are  deter¬ 
mined  by  a  complex  interplay  between  a  number  of  different  mechanisms.  At  short 
residence  times,  the  carbon  fractions  in  the  form  of  carbon  monoxide,  CHm,  and 
C2Hn  species  can  depend  strongly  on  the  dissociation  pathways  of  the  carbon-bearing 
source  gas.  As  an  example,  we  consider  the  chemistry  of  ethanol.  Hydrogen  atoms 
attack  ethanol  via  four  important  reactions.  Two  of  these, 

C2H5OH  +  H  =  sC2H4OH  +  H2,  (16) 

C2H5OH  +  H  =  C2H50  +  H2,  (17) 

lead  directly  to  the  formation  of  carbon  monoxide  and  methyl  radicals  via  the  follow¬ 


ing  paths 

sC2H4OH  -»  CH3CHO  ->  HCO  -4  CO,  (18) 

sC2H4OH  -4  CH3CHO  -4  CH3CO  -4  CO,  (19) 

C2H50  ->  CH20  -4  HCO  ->  CO.  (20) 

A  second  pair  of  reactions, 

C2H5OH  +  H  =  PC2H4OH  +  H2,  (21) 

C2H5OH  +  H  =  C2H5  +  H20,  (22) 


result  in  the  formation  of  C2  species,  either  directly  or  through  a  pathway  involving 
the  unimolecular  dissociation 

pC2H4OH  =  C2H4  +  OH.  (23) 
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Since  all  four  of  the  initial  reactions  in  the  ethanol  decomposition  involve  the  same 
species,  the  branching  ratio  between  the  formation  of  Ci  and  C2  species  depends 
solely  on  the  relative  Arrhenius  rates.  The  reaction  rates  and  branching  ratio  are 
shown  as  a  function  of  temperature  in  Fig.  9.  As  the  temperature  is  raised  above 
2000  K,  the  formation  of  C2  species  becomes  favored,  but  it  should  be  noted  that  over 
the  range  1000-4000  K,  the  two  pathways  are  comparable  with  the  branching  ratio 
only  varying  from  0.7  to  1.3. 

Although  the  dissociation  pathways  of  the  source  gases  can  have  a  strong  effect 
on  the  carbon  and  hydrogen  distributions,  the  generic  elementary  reactions  common 
to  all  carbon-hydrogen-oxygen  mixtures  are  usually  more  important.  Under  the  con¬ 
ditions  of  interest  in  this  study,  the  dominant  pathways  for  the  creation  of  C2  species 
from  Ci  species  are  the  reactions  of  the  forms 


CHn  +  CHm  =  C2Hn+m_i  +  H,  (24) 

CHn  +  CHm  =  C2Hn+m_2  4-  H2.  (25) 

Since  these  are  simple  two-body  reactions,  they  can  proceed  much  more  readily  at 
low  gas  densities  than  the  pressure-dependent,  three-body  reaction 

CH3  +  CH3(+M)  =  C2H6(+M),  (26) 

which  is  important  for  methane  under  typical  combustion  conditions.  Reactions  of 
the  more  general  forms 

CnHm  +  CpHq  =  Cn+pHm+q_i  +  H,  (27) 

C„Hm  +  CpHq  =  Cn+pHm+q_2  +  H2,  (28) 
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are  responsible  for  the  creation  of  heavier  hydrocarbon  species. 

While  only  the  forward  directions  of  reactions  (25)  and  (28)  are  generally  im¬ 
portant,  the  reverse  reactions  (-24)  and  (-27)  play  an  important  role  in  the  pyrolysis 
of  heavier  hydrocarbons.  The  balance  between  heavy  and  light  hydrocarbon  species 
depends  strongly  on  both  the  reaction  rates  and  the  hydrogen  atom  concentration. 
The  fraction  of  hydrogen  in  the  form  of  hydrogen  atoms  in  turn  depends  on  the  rates 
of  the  dissociation  reactions 


h2  +  m  =  h  +  h  +  m, 

(29) 

RH  +  M  =  R  +  H  +  M, 

(30) 

and  the  hydrogen  shuffling  reactions 

H2  +  R  =  RH  +  H, 

(31) 

where  R  is  a  generic  hydrogen-carbon-oxygen  species.  The  hydrogen  atom-carbon 
radical  ratio  is  also  affected  by  the  creation  of  carbon  monoxide.  Carbon  monoxide 
can  be  formed  either  directly  from  the  dissociation  of  the  source  gases,  as  described 
above  for  ethanol,  or  through  the  reactions  involving  oxygen  atoms  or  molecules, 
hydroxyl  radicals,  and  hydrocarbon  species. 

The  carbon  and  hydrogen  profiles  in  Figs.  2-7  clearly  show  the  interplay  between 
these  various  mechanisms.  As  an  example,  we  consider  the  CO-rich  methane-oxygen 
mixture.  For  this  source  gas  recipe,  methane  is  the  only  source  of  hydrogen.  As 
a  result,  in  the  limit  of  very  high  flow  rates,  Ci  radicals  and  hydrogen  atoms  are 
formed  at  almost  the  same  rate.  Decreasing  the  flow  rate  allows  the  fast  reactions 
(24,25,27,28)  to  occur,  leading  to  the  formation  of  heavier  hydrocarbon  species  from 
the  Ci  radicals.  As  the  flow  rate  is  further  decreased,  slower  reactions  which  result 
in  the  formation  of  hydrogen  atoms  and  carbon  monoxide  become  more  important 
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and  the  gas  mixture  asymptotically  approaches  the  thermodynamic  equilibrium  com¬ 
position.  The  minimum  in  /( Cirad)  at  approximately  m  =  0.01  g/s  is  due  to  the 
competition  between  the  forward  reactions  (24,25,27,28)  which  consume  Ci  radicals 
at  high  flow  rates  when  /(H)  is  low  and  the  reverse  reactions  (-24,-27)  which  produce 
Ci  radicals  at  low  flow  rates  when  /(H)  is  high.  The  behavior  of  Dmax  in  the  limit 
777  — ►  oo  follows  from  the  fact  that  the  formation  of  methyl  radicals  at  short  residence 
times  is  due  almost  entirely  to  the  dissociation  reaction 

CH4  +  M  =  CH3  +  H  +  M.  (32) 

In  the  high-flow-rate  limit,  /(Cirad)  becomes  directly  proportional  to  residence  time, 
or  inversely  proportional  to  flow  rate,  resulting  in  Dmax  independent  of  flow  rate.  Of 
course,  since  the  ratio  of  hydrogen  atoms  to  methyl  radicals  approaches  unity,  the 
growth  of  high-quality  diamond  would  not  be  expected. 

While  many  of  the  features  described  above  for  the  methane-oxygen  mixture  are 
common  to  the  other  gas  recipes  studied,  there  are  important  qualitative  differences. 
For  the  gas  mixtures  utilizing  ethanol  or  acetylene  as  the  carbon  source,  Dmax  exhibits 
a  global  maximum  as  a  function  of  flow  rate  and  decreases  in  the  limit  m  —*■  oo.  This  is 
because  the  initial  step  in  the  dissociation  of  the  source  gases  does  not  directly  yield 
Ci  radicals,  but  rather  intermediate  species  which  must  undergo  further  reactions 
before  producing  growth  species. 

The  analysis  up  to  this  point  has  been  based  on  the  assumption  that  the  disso¬ 
ciation  of  hydrogen  is  due  entirely  to  gas-phase  reactions  involving  neutral  species. 
In  reality,  the  situation  is  typically  more  complex  since  H2  molecules  can  dissociate 
on  the  filament  surfaces  in  hot-filament  reactors  or  through  electron-impact  reactions 
in  plasma-enhanced  CVD.  To  account  for  these  additional  mechanisms  of  hydrogen 
dissociation,  we  have  run  a  set  of  simulations  in  which  the  rate  k^iss  °f  the  dissoci- 
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ation  reaction  (29)  has  been  increased  by  factors  of  10,  100,  and  1000.  The  results 
for  T  and  Dmax  obtained  using  the  ethanol-water-hydrogen  and  ethanol-hydrogen 
mixtures  are  displayed  in  Figs.  10  and  11.  Both  $  and  Dmax  generally  increase  for  a 
given  flow  rate  as  is  increased,  with  gains  of  up  to  an  order  of  magnitude  at  the 
largest  values  of  kdlss.  This  has  the  effect  of  pushing  the  maximum  flow  rate  at  which 
diamond  growth  is  expected  to  occur  (T  >  104),  and  hence  the  maximum  diamond 
growth  rate,  to  higher  values.  The  peak  in  T  at  flow  rates  of  approximately  0.01 
g/s  is  seen  to  diminish  or  even  vanish  as  kdiss  is  increased.  This  is  due  to  the  fact 
that  the  increased  hydrogen  dissociation  drives  the  system  towards  thermodynamic 
equilibrium  more  quickly  and  reduces  or  eliminates  the  dip  in  /(Cirad). 

V.  SUMMARY 

The  Bachmann  diagram  has  proven  to  be  a  very  useful  construct  in  predicting 
which  source-gas  recipes  are  likely  candidates  for  diamond  growth,  but  its  utility 
may  be  limited  because  it  does  not  account  for  important  kinetic  effects  in  the  gas- 
phase  chemistry.  The  fact  that  the  diamond  growth  regime  is  determined  solely  by 
the  overall  C-H-0  stoichiometry  suggests  that  the  source  gases  have  been  allowed 
to  react  for  a  relatively  long  time  and  establish  a  quasi-thermodynamic  equilibrium 
before  reaching  the  growth  surface.  In  most  diamond  CVD  reactors  this  is  the  case. 
In  low-pressure,  high-flow-rate  reactors  where  the  chemical  reaction  time  scales  are 
long  or  comparable  to  the  times  required  for  diffusion  to  the  growth  surface  or  flow 
through  the  reactor  to  occur,  however,  this  assumption  may  no  longer  be  true  and 
kinetic  effects  must  be  taken  into  account. 

In  this  paper,  we  have  studied  the  effect  of  varying  the  source-gas  flow  rates  on 
both  the  diamond  quality  and  the  theoretical  maximum  growth  rate.  We  have  found 
that  increasing  the  gas  flow  rate  generally  increases  the  growth  rate,  but  often  at  the 
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expense  of  material  quality.  At  high  flow  rates,  there  can  be  both  an  unfavorable  ratio 
of  hydrogen  atoms  to  Ci  radicals,  which  we  have  quantified  in  terms  of  Goodwin  s 
quality  parameter  [6],  and  large  concentrations  of  C3  and  heavier  hydrocarbon  species, 
resulting  in  degradation  of  the  diamond  film.  For  gases  utilizing  a  C2  carbon  source, 
the  maximum  growth  rate  has  been  shown  to  decrease  at  high  flow  rates  since  the 
residence  time  in  the  reactor  is  too  short  to  allow  for  the  formation  of  potential  growth 
species.  Important  quantitative  differences  exist  between  the  CO-rich  and  CO-lean 
mixtures.  Because  the  CO-rich  systems  contain  large  amounts  of  carbon  not  initially 
bound  in  the  form  of  carbon  monoxide,  unfavorably  high  concentrations  of  C3  and 
heavier  hydrocarbon  species  may  be  built  up  at  high  flow  rates. 

Of  the  various  recipes  studied,  the  CO-lean,  acetylene-oxygen-hydrogen  mixture 
appears  to  be  the  most  robust.  For  all  flow  rates,  the  atomic  hydrogen-to-carbon 
radical  ratio  is  predicted  to  be  sufficiently  large  so  as  to  yield  high-quality  diamond 
growth.  Unfortunately,  this  result  is  apparently  of  theoretical  interest  only,  since  the 
explosive  hazards  associated  with  using  acetylene-hydrogen  mixtures  in  a  microwave- 
plasma  or  hot-filament  reactor  far  outway  the  potential  benefits. 

Obviously,  the  perfectly-stirred  reactor  calculations  ignore  many  important  phys¬ 
ical  processes  that  occur  within  a  real  diamond  CVD  reactor.  Surface  reactions  will 
have  the  effect  of  both  depleting  the  carbon  radical  population  and  converting  atomic 
to  molecular  hydrogen.  Numerical  simulations  taking  into  account  plasma  chemistry 
indicate  that  the  electron  impact  dissociation  of  carbon  monoxide  may  be  impor¬ 
tant,  resulting  in  a  considerably  larger  flux  of  growth  species  to  the  surface  [34].  The 
counter-diffusion  of  molecular  hydrogen  and  atomic  hydrogen  through  the  boundary 
layer  over  the  growth  surface  has  also  been  ignored.  The  study  of  these  effects  is  left 
to  future  work.  Nonetheless,  the  general  results  derived  from  this  study  regarding  the 
dependence  of  diamond  quality,  maximum  growth  rate,  and  gas-phase  composition 
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on  the  source  gases  and  flow  rates  are  expected  to  be  approximately  valid. 
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Fig.  1.  C-H-0  Ternary  diagram.  Dashed  lines,  adapted  from  ref.  [1],  indicate  approximate 
locations  of  boundaries  separating  "diamond  growth"  domain  from  "no  growth"  and  "non¬ 
diamond  growth"  domains.  Dotted  line  is  the  CO-line  and  represents  stoichiometries  with  equal 
mole  fractions  of  carbon  and  oxygen.  Open  circles  represent  locations  of  pure  compounds  and 
filled  circles  indicate  CO-rich  (mixl)  and  CO-lean  (mix2)  stoichiometries  investigated  in  this 
work. 
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Fig.  2.  Results  of  PSR  calculations  for  ethanol-water-hydrogen  21:20:1.  (a)  Quality  factor,/ 
(Qrad),  /(H),  and  maximum  growth  rate  (dotted  line,  plotted  against  right  axis)  versus  flow 
rate,  (b)  carbon  and  hydrogen  distributions  versus  flow  rate. 
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Fig.  2.  (continues)  Results  of  PSR  calculations  for  ethanol-water-hydrogen  21:20:1.  (a)  Quality 
factor, /(Cjrad), /(H),  and  maximum  growth  rate  (dotted  line,  plotted  against  right  axis)  versus 
flow  rate,  (b)  carbon  and  hydrogen  distributions  versus  flow  rate. 
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Fig.  3.  Results  of  PSR  calculations  for  methane-oxygen  84:41.  (a)  Quality  factor,  /(Qrad),/ 
(H),  and  maximum  growth  rate  (dotted  line,  plotted  against  right  axis)  versus  flow  rate,  (b) 
carbon  and  hydrogen  distributions  versus  flow  rate. 


Fig.  3.  (continues)  Results  of  PSR  calculations  for  methane-oxygen  84:41.  (a)  Quality  factor, 
/(Qrad),/^),  and  maximum  growth  rate  (dotted  line,  plotted  against  right  axis)  versus  flow 
rate,  (b)  carbon  and  hydrogen  distributions  versus  flow  rate. 
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Fig.  4.  Results  of  PSR  calculations  for  acetylene-oxygen-hydrogen  42:41:126.  (a)  Quality  factor, 
/(Qrad^/CH),  and  maximum  growth  rate  (dotted  line,  plotted  against  right  axis)  versus  flow 
rate,  (b)  carbon  and  hydrogen  distributions  versus  flow  rate. 


Fig.  4.  (continues)  Results  of  PSR  calculations  for  acetylene-oxygen-hydrogen  42:41:126.  (a) 
Quality  factor, /(Qrad), /(H),  and  maximum  growth  rate  (dotted  line,  plotted  against  right  axis) 
versus  flow  rate,  (b)  carbon  and  hydrogen  distributions  versus  flow  rate. 
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Fig.  5.  Results  of  PSR  calculations  for  ethanol-hydrogen  1:81.  (a)  Quality  factor,  / (Qrad),  / 
(H),  and  maximum  growth  rate  (dotted  line,  plotted  against  right  axis)  versus  flow  rate,  (b) 
carbon  and  hydrogen  distributions  versus  flow  rate. 
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Fig.  5.  (continues)  Results  of  PSR  calculations  for  ethanol-hydrogen  1:81.  (a)  Quality  factor, 
/ (Cirad),/(H),  and  maximum  growth  rate  (dotted  line,  plotted  against  right  axis)  versus  flow 
rate,  (b)  carbon  and  hydrogen  distributions  versus  flow  rate. 
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Fig.  6.  Results  of  PSR  calculations  for  methane-oxygen-hydrogen  4:1:160.  (a)  Quality  factor, 
/(C,rad),/(H),  and  maximum  growth  rate  (dotted  line,  plotted  against  right  axis)  versus  flow 
rate,  (b)  carbon  and  hydrogen  distributions  versus  flow  rate. 


Fig.  6.  (continues)  Results  of  PSR  calculations  for  methane-oxygen-hydrogen  4:1:160.  (a) 
Quality  factor, /(Qrad),  /(H),  and  maximum  growth  rate  (dotted  line,  plotted  against  right  axis) 
versus  flow  rate,  (b)  carbon  and  hydrogen  distributions  versus  flow  rate. 
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Fig.  7.  Results  of  PSR  calculations  for  acetylene-oxygen-hydrogen  2:1:166.  (a)  Quality  factor, 
/(Ct rad), /(H),  and  maximum  growth  rate  (dotted  line,  plotted  against  right  axis)  versus  flow 
rate,  (b)  carbon  and  hydrogen  distributions  versus  flow  rate. 
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Fig.  7.  (continues)  Results  of  PSR  calculations  for  acetylene-oxygen-hydrogen  2:1:166.  (a) 
Quality  factor, /(Qrad), /(H),  and  maximum  growth  rate  (dotted  line,  plotted  against  right  axis) 
versus  flow  rate,  (b)  carbon  and  hydrogen  distributions  versus  flow  rate. 


35 


1  o'5  1  o'4  1  o'3  1  0‘2  1  o’1 ,  1,0°  1  o’  1  o2  1  o3 
flow  rate  (g/s) 


Fig.  8.  Quality  factors  for  the  six  different  source-gas  mixtures  versus  flow  rate 
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Fig.  9.  Arrhenius  rates  for  hydrogen  atom-ethanol  reactions  yielding  products  (a)  pC2H4OH,  (b) 
sC2H4OH,  (c)  C2H50  or  C2H5,  and  branching  ratio  (plotted  against  right  axis)  between  reactions 
leading  to  Q  and  C2  hydrocarbon  formation  versus  temperature. 
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Fig.  10.  (continues)  Results  of  PSR  calculations  for  ethanol- water-hydrogen  21:20:1  at  various 
H2  dissociation  rates,  (a)  Quality  factors  and  (b)  maximum  growth  rates. 


Fig.  11.  Results  of  PSR  calculations  for  ethanol-hydrogen  1:81  at  various  H2  dissociation  rates, 
(a)  Quality  factors  and  (b)  maximum  growth  rates. 
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Fig.  11.  (continues)  Results  of  PSR  calculations  for  ethanol-hydrogen  1:81  at  various  H2 
dissociation  rates,  (a)  Quality  factors  and  (b)  maximum  growth  rates. 
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